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The superwind mass-loss rate of the metal-poor carbon 
star LI-LMC 1813 in the LMC cluster KMHK 1603* 
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ABSTRACT 

LI-LMC 1813 is a dust-enshrouded Asymptotic Giant Branch (AGB) star, located in 
the small open cluster KMHK 1603 near the rim of the Large Magellanic Cloud (LMC). 
Optical and infrared photometry between 0.5 and 60 fixa is obtained to constrain the 
spectral energy distribution of LI-LMC 1813. Near-infrared spectra unambiguously 
show it to be a carbon star. Modelling with the radiation transfer code Dusty yields 
accurate values for the bolometric luminosity, L — 1.5 x 10 4 L Q , and mass-loss rate, 
M = 3.7(±1.2) x 10~ 5 M© yr _1 . On the basis of colour-magnitude diagrams, the age 
of the cluster KMHK 1603 is estimated to be t = 0.9-1.0 Gyr, which implies a Zero- 
Age Main Sequence mass for LI-LMC 1813 of A/ Z ams = 2.2 ± 0.1 M . This makes 
LI-LMC 1813 arguably the object with the most accurately and reliably determined 
(circum)stellar parameters amongst all carbon stars in the superwind phase. 
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1 INTRODUCTION 

In the final stages of their evolution, stars of main-sequence 
masses Mzams ~ 0.8-8 Mq ascend the Asymptotic Giant 
Branch (AGB). The end-product of their evolution is a 
White Dwarf, with a typical mass of Mwd ~0.5 Mq. It 
is thus clear that near the tip of the AGB, stars must be 
losing most of their mass. The rate at which this occurs is 
estimated to be M ~ 10 -6 to 10 -4 Mq yr -1 , much faster 
than the nuclear burning mass consumption rate which nor- 
mally determines the rate of stellar evolution (van Loon et 
al. 1999b). Hence the mass loss is important for the sub- 
sequent evolution of the star, as well as for the chemical 
enrichment of the inter-stellar medium. Yet our knowledge 
and understanding of the mass-loss mechanism and the com- 
position of the mass lost is incomplete. 

Amongst the key factors in supporting a dense outflow 
from an AGB star are the low effective temperature, which 
allows for the formation of circum-stellar dust, and the high 
luminosity, which, through radiation pressure on the dust 
grains and collissional coupling between the dust and the 
gas, drives a stellar wind. Because the dust only forms at 
several stellar radii from the star, there has to be an addi- 
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tional mechanism to bring the gas into the dust formation 
region — possibly provided by strong radial pulsations (e.g. 
Bowen & Willson 1991). 

Dynamical models for the circumstellar envelopes of 
AGB stars suggest that the mass-loss rate depends strongly 
on the stellar effective temperature, about MtxT^ (Arndt, 
Fleischer & Sedlmayr 1997) . Although this seems to be con- 
firmed for the onset of the dusty mass-loss epoch (Alard et 
al. 2001), it is not clear whether such strong temperature 
dependence is also valid during the so-called "superwind" 
phase where the AGB star is obscured by its optically thick 
circum-stellar envelope (van Loon et al. 1999b). 

A dramatic change in the chemical composition of the 
mass lost occurs when an AGB star becomes a carbon star 
(Iben & Renzini 1983): during the temporal ignition of a 
helium-burning shell, situated in between the degenerate 
carbon/oxygen core and the hydrogen- burning shell, the 
convective mantle may reach the nuclear burning layer and 
transport nuclear burning products to the stellar surface (3 rd 
dredge- up). These thermal pulses (TP) occur on timescales 
of tTP ~ 10 4 -10 5 yr, and last only briefly. But over the 
course of TP- AGB evolution, the surface abundance of car- 
bon may exceed that of oxygen, forming a carbon star. Due 
to the strong binding of the CO molecule, only the remain- 
ing oxygen or carbon is used to synthesize other molecules 
that nucleate to form dust grains in the circum-stellar enve- 
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lope. Hence the circum-stellar chemistry of carbon stars is 
fundamentally different from that of oxygen-rich AGB stars. 

There is a threshold mass above which 3 rd dredge- 
up is strong enough to make a carbon star, roughly at 
Mzams > 1-2 Mq and mildly dependent on metallicity 
(Marigo, Girardi & Bressan 1999). However, in massive AGB 
stars, about Mzams > 4 Mq, nuclear burning occurs within 
the bottom layers of the convective mantle, converting the 
carbon into nitrogen: Hot Bottom Burning (HBB). At the 
very end of AGB evolution, the diminished mantle mass as 
a result of prolonged mass loss causes HBB to cease, and 
a final thermal pulse may turn the star into a carbon star 
after all (Frost et al. 1998). These competing processes are 
likely to depend on the initial metallicity of the star. 

Optical studies of the AGB population in the Large 
Magellanic Cloud (LMC) provide support for the 3 rd dredge- 
up and HBB to govern the formation of carbon stars (e.g. 
Groenewegen & de Jong 1993; Marigo et al. 1999). Recent 
infrared (IR) studies could assess the effects of mass loss and 
possible observational bias in optical studies. These studies 
confirm the delayed formation of massive carbon stars after 
cessation of HBB, but also suggest that amongst AGB stars 
in the superwind phase oxygen-rich stars exist with masses 
below the threshold for HBB (Zijlstra et al. 1996; van Loon 
et al. 1997, 1998, 1999b; van Loon, Zijlstra & Groenewegen 
1999a). 

Much of the uncertainty with respect to the progenitor 
mass and initial metallicity of an individual AGB star can be 
alleviated if the star is a member of a co-eval stellar cluster. 
Optical studies provide a great number of such examples 
in LMC and SMC clusters, confirming the mass-range for 
optically bright carbon stars between Mzams, min ~ 1—1.5 
and -Mzams, max ~3-5 M Q (Frogel, Mould & Blanco 1990; 
Marigo, Girardi & Chiosi 1996; van Loon 2002). 

Van Loon et al. (1998) noted that the oxygen-rich dust- 
enshrouded OH maser star IRAS 05298—6957 is located in 
the LMC cluster HS 327, and their discovery of an opti- 
cally bright carbon star in the same cluster directly provided 
a value for the threshold mass above which HBB occurs: 
Mzams = 4 M (van Loon et al. 2001). 

Tanabe et al. (1997) found one dust-enshrouded carbon 
star in each of the LMC clusters NGC 1783 and NGC 1978, 
and another one in the SMC cluster NGC 419. They all have 
moderate luminosities, 5000 < L < 8000 Lq, and progeni- 
tor masses, Mzams ~ 1.5-1.6 Mq. Their mass-loss rates are 
M < 10 Mq yr~ , but the relatively long pulsation pe- 
riods of P ~ 500 d indicate that they have already shed a 
significant fraction of their mass (Nishida et al. 2000). 

The IRAS point source LI-LMC 1813 (IRAS 
06025—6712) was identified as a candidate dust-enshrouded 
AGB star in the LMC by Loup et al. (1997). Its suspected 
nature was confirmed after the identification of the ex- 
tremely red counterpart (van Loon et al. 1997), and it was 
later found to be coincident with an LMC cluster (van Loon 
1999), KMHK 1603 (Kontizas et al. 1990). Kontizas et al. 
estimate that the number density of LMC star clusters in 
the direction of LI-LMC 1813 is ~20 per square degree, 
with each cluster covering, on average, less than 1 square 
arcminute. Hence, in the direction of LI-LMC 1813, less 
than 0.6% of sky is covered by an LMC cluster. There 
are of order 10 2 IRAS-selected AGB stars in the LMC 
(Loup et al. 1997; van Loon et al. 1997), most of which are 
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Figure 1. Image of KMHK 1603 and surrounding field, taken 
with the 0.9m Dutch telescope through the i-band filter. The 
brightest star near the centre of the field is KMHK 1603-1, whilst 
LI-LMC 1813 is invisible at this wavelength. 



concentrated towards the LMC bar. One does not therefore 
expect any chance coincidence of an IRAS-selected AGB 
star within the < 0.6% of sky covered by an LMC cluster 
in this part of the LMC, rendering LI-LMC 1813 likely to 
be a physical member of the LMC cluster KMHK 1603. 

We present here the results of extensive optical and IR 
imaging and spectroscopic follow-up observations of LI-LMC 
1813 and the cluster KMHK 1603. We identify the dust- 
enshrouded AGB star as a carbon star, and derive accurate 
estimates for its progenitor mass, Mzams = 2.2 Mq, and 
mass-loss rate, M = 3.7 x 10~ 5 Mq yr _1 , which we confront 
with theoretical models. Remarkably, the mass- loss rate is 
very high for this metal-poor (~ | solar) object. 



2 THE BRIGHTNESS OF LI-LMC 1813 

We have collected optical (Gunn g, r & i), near-IR (J, K & 
L') and mid-IR (N) imaging data, which we combine with 
archival IR data (2MASS, Denis, MSX and IRAS) to de- 
scribe the spectral energy distribution (SED) of LI-LMC 
1813. These data are described in full in Appendix Al. 



2.1 Images 

The moderately crowded optical images (Fig. 1) feature the 
sparse cluster KMHK 1603, dominated by the bright star 
KMHK 1603-1. At these wavelengths, LI-LMC 1813 is in- 
visible. Undetected in the J-band, the brightness of LI-LMC 
1813 in the K-band rivals that of KMHK 1603-1 (Fig. 2). 
These are the only two stars visible at A ~ 3 /jm, with LI- 
LMC 1813 by then outshining KMHK 1603-1. At even longer 
wavelengths, only LI-LMC 1813 is detected. 
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Table 1. Photometry of LI-LMC 1813 and the optically bright star KMHK 1603-1: magnitudes for the optical and near-IR, and flux 
densities (in Jy) for the mid-IR (A > 6 (im). For Denis, MSX and IRAS no exact epoch is available. 
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Figure 2. Near-IR image of KMHK 1603, taken with IRAC-2. 
The brightest star in the field, near the edge, is KMHK 1603-1. 
The very red IR object LI-LMC 1813 is almost equally bright in 
the K-band (left; arrow) but invisible in the J-band (right). 
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Figure 3. K-band (solid disks) and L'-band (solid squares) pho- 
tometry of LI-LMC 1813, with respect to their mean magnitudes, 
as a function of Julian Date. The data can be matched by a si- 
nusoidal curve with period P = 680 d. The open symbols are the 
data repeated each 680 d. 



2.2 Photometry 

The complete set of optical and infrared photometry in 
the range 0.5-60 ^m for LI-LMC 1813 is listed in Table 1. 
Also included in the table is photometry for KMHK 1603-1, 
whose IRAC-2, 2MASS and Denis magnitudes, which have 
been transformed onto the SAAO system (Carter 1990; Car- 
penter 2001), are in excellent agreement with each other. 

2.3 Variability 

A third epoch of K-band photometry is available from the 
acquisition image for the H+K-band spectroscopy in July 
2002 (Appendix A2). The K-band magnitude of LI-LMC 
1813 is obtained by assuming that the K-band magnitude of 
KMHK 1603-1 remained unchanged at K = 11.80 mag. 



A second epoch of L'-band photometry is obtained from 
the acquisition image taken for the L-band spectroscopy in 
December 2001 (Appendix A2). The L'-band magnitude of 
KMHK 1603-1 is assumed to have remained unchanged at 
L' = 11.75 mag (if the zero point were the same as that 
in October 2000, the L'-band magnitude of KMHK 1603-1 
would be L' = 11.80 ± 0.30 mag). Hence the L'-band mag- 
nitude of LI-LMC 1813 is V = 9.27 ±0.25 mag, in excellent 
agreement with the flux-calibrated spectrum. 

A sparse lightcurve is constructed from the differences 
of the three K-band and two L'-band magnitudes with their 
mean (Fig. 3): the variability of (obscured) AGB stars in 
the K and L'-bands are in phase with each other, and have 
similar amplitudes (see Whitelock et al. 2003). The data 
are consistent with a pulsation period of P ~ 680 d. If con- 
firmed, this would be rather long but not unusual for highly 
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evolved AGB stars (Wood 1998) . The large L'-band ampli- 
tude of AL' ~1.6 mag might be due in part to temporal 
variations in the amount and temperature of the warm dust 
contributing to the IR excess at A =3-4 /im. 



2.4 Interstellar extinction 

The interstellar extinction towards the LMC is generally 
low (but see van Loon et al. 1997 for severely reddened 
background galaxies) . From the reddening maps provided by 
Burstein & Heiles (1982), based on galaxy counts, the red- 
dening towards LI-LMC 1813 is estimated to be E(B — V) = 
0.06 mag. The reddening maps derived from the far-IR dust 
emission mapped by COBE/DIRBE (Schlegel, Finkbeiner 
& Davis 1998) yield a very consistent value of E(B — V) = 
0.056 mag, which we adopt here. The spatial resolution of 
these maps is insufficient to detect differences in redden- 
ing over a few arcminutes, and therefore uniform reddening 
is assumed for the region within a radius of ~ 2' around 
KMHK 1603. This converts to an extinction of A\ = 0.20, 
0.15 and 0.11 in the g, r and i-bands, and A\ — 0.05, 0.03, 
0.02 and 0.01 mag in the J, H, K and L'-bands, respectively 
(Mathis 1990). The (uncorrected) photometry of Table 1 is 
corrected for this extinction before further analysis. 



2.5 The distance to LI-LMC 1813 

For the distance modulus of the LMC we adopt (m — M)o = 
18.55±0.17 mag as derived by Groenewegen (2000) from the 
Cepheid period-luminosity relationship. However, the LMC 
is inclined, causing the distance modulus to vary across the 
face of the LMC. This is especially important for LI-LMC 
1813, which is situated near the rim of the galaxy. Using the 
values for the centre, inclination and position angle of the 
nodes from van der Marel & Cioni (2001), the differential 
distance modulus with respect to the centre of the LMC is 
A(m — M) = —0.08 mag. Within an uncertainty of ~ 0.2 
mag, the distance modulus to LI-LMC 1813 then becomes 
(m - M) = 18.47 mag. 



3 THE PROGENITOR MASS OF LI-LMC 1813 

The progenitor mass of LI-LMC 1813 may be estimated 
by isochrone-fitting to the colour-magnitude diagram of the 
stars in the cluster KMHK 1603. 



3.1 Multi-object photometry 

Multi-object photometry was performed on the optical im- 
ages using an implementation of DAOPHOT (version n) 
and ALLSTAR (Stetson 1987) within ESO-MIDAS. The 
reader is referred to van Loon et al. (2001) for a detailed 
description of the Point Spread Function (PSF), multi- 
object photometry and photometric calibration. As a result, 
~ 7 x 10 3 sources were detected in each of the gri-bands in 
the 3.7 x 10~ 3 square degree area around KMHK 1603. The 
completeness drops rapidly beyond ~ 23 mag. 

The point-sources were cross-correlated between the dif- 
ferent filters after determining and applying a geometric 
transformation (rotation and linear translation), using an 
iterative scheme with a growing search radius and rejecting 
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Figure 4. Stellar surface density in the g-band. Contour levels 
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extended sources on the basis of the sharpness parameter re- 
turned by ALLSTAR. The photometry was then corrected 
for interstellar extinction as derived in section 2.4. 

Because there are only ~ 10 2 objects visible in the near- 
IR images, and hence these images are not at all crowded, 
the near-IR photometry was obtained manually by aper- 
ture photometry. The photometry was transformed onto the 
SAAO system (see Appendix A1.2) and then corrected for 
interstellar extinction (section 2.4). 

The surface density of stars in the g-band image was 
measured by counting stars with g < 22 mag and sharpness 
s G [—0.5,0.5], within a 28" x 28" box which was succes- 
sively moved across the image at steps of ~ 9". The result- 
ing contour map (Fig. 4) clearly shows the cluster KMHK 
1603 as a region with a higher density of stars as compared 
to the surrounding field. The cluster appears to be some- 
what elliptical, but for practical purposes we assume it to 
be spherical. The centre of the cluster is determined to be 
located at (J2000) RA= 6 h 02 m 29 s and Dec= -67°12'54". 

The surface density of stars, as measured in 6" annuli 
centred on the cluster, is given as a function of the radial 
distance in Fig. 5. The surface density in the core of the 
cluster is three times that of the field, and it drops to the 
field level at a radius of 40". With a (projected) distance 
to the cluster centre of 15", LI-LMC 1813 is located well 
within the cluster boundaries. To minimise contamination 
by field stars, the colour-magnitude diagrams of the cluster 
were constructed for stars within a smaller radius of 30". 



3.2 The age and metallicity of KMHK 1603 

The interstellar extinction-corrected io versus (g—i)o colour- 
magnitude diagram is presented for the cluster KMHK 1603 
(Fig. 6, left) and the surrounding field (Fig. 6, right). Also 
plotted are three isochrones from Bertelli et al. (1994), for a 
typical LMC metallicity of [M/H]= -0.4 and ages of 10 Gyr, 
1 Gyr and 100 Myr, after transformation onto the photomet- 
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Figure 6. Colour-magnitude diagram of iq versus (g—i)o for the LMC cluster KMHK 1603 (left) and surrounding field (right). Isochroncs 
(Bertelli et al. 1994) are plotted for [M/H]= -0.4, and ages of 10 8 , 10 9 and 10 10 yr. LI-LMC 1813 is not detected at these wavelengths. 
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Figure 5. Radial stellar surface density profile in the g-band 
(solid: smoothed; dashed: unsmoothed), with respect to the centre 
of the cluster KMHK 1603. 
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Figure 7. Close-up of the Main-Sequence turnoff and red clump 
regions of the colour-magnitude diagram of io versus (g — i)o for 
the LMC cluster KMHK 1603. Isochrones (Bertelli et al. 1994) 
arc plotted for several combinations of age and metallicity. 



ric systems as described in detail in van Loon et al. (2003). 
These isochrones do not include the effects of circumstellar 
extinction. 

Both the cluster and field show a well-populated Main 
Sequence up to an age of t ~ 10 9 yr, and a red clump of 
core He-burning giants, at (g — i)o ~ 1 and io ~ 18 mag, 
consistent with such intermediate or older ages. Although 



there are a number of field Main Sequence stars brighter 
than the 10 9 yr Main Sequence turnoff, and corresponding 
RGB or AGB stars with io ~ 16 mag, there is a clear absence 
of stars younger than a few 10 8 yr. 

A close-up of the optical colour-magnitude diagram for 
the cluster stars (Fig. 7) reveals two populations: one of t ~ 1 
Gyr and an older one of t ~ 2 Gyr. The younger population is 
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tions (Fig. 8), but it is clear that the cluster population is 
not consistent with an age of much less than t ~ 1 Gyr (note 
that the 2MASS data sample a much larger area than the 
IRAC-2 data). The red clump is clearly visible in the cluster, 
at (J — K)o ~ 0.5 and Kq ~ 17 mag. 

The brightness of the red clump stars can be used as 
an age indicator. In the optical, the red clump in the cluster 
KMHK 1603 is observed at rru = 18.10 ± 0.20 mag (Fig. 
7) or, assuming that mi~mi, at Mi = —0.37 ± 0.20 mag. 
This suggests an age of either t = 0.82l° ^ Gyr or t >1.22 
Gyr for [M/H]= -0.4, and t = 0.96lg° 16 Gyr for [M/H] = 
-0.7 (Girardi & Salaris 2001, their Table 1). In the ncar- 
IR, the red clump is observed at ttik = 16.90 ± 0.20 mag, 
or Mk = — 1.57 ± 0.20 mag. This suggests an age of either 
t = 0.72l2:i| Gyr or t >1.32 Gyr for [M/H]= -0.4, and 
either t = 6.78lg;^ Gyr or t >1.39 Gyr for [M/H]= -0.7 
(Salaris & Girardi 2002, their Table 1). The ambiguity in 
the age solution is the reason for the near-coincidence of the 
red clump for the cluster and field stars, if the cluster and 
field are dominated by populations of t ~ 1 Gyr and t > 1.4 
Gyr, respectively. 



Figure 8. Colour-magnitude diagram of Kq versus (J — K)o 
for the LMC cluster KMHK 1603 (solid) and surrounding field 
(open), both from the IRAC-2 imaging as well as from the 2MASS 
database. Isochrones (Bertelli et al. 1994) arc plotted for several 
combinations of age and metallicity. The AGB star LI-LMC 1813 
is extremely red due to circumstellar extinction. 



3.3 Conclusion: progenitor mass of LI-LMC 1813 

We conclude that LI-LMC 1813 probably has an age of 
t = 1 Gyr (or slightly younger) and an initial metallicity 
of [M/H]= —0.7 (or slightly higher), and hence a progenitor 
mass of Mzams = 2.2 ± 0.1 M (Bertelli et al. 1994). 



traced by a Main Sequence turnoff (MS-TO) at (g — i)o ~ 0.2 
and io ~ 19.0 mag, clearly separated from the older popula- 
tion with a MS-TO at (g-i)o ~ 0.4 and i ~ 19.8 mag. Both 
populations show a red clump at roughly the same location 
in the colour-magnitude diagram. The younger population 
stands out prominently amongst the stars in the direction of 
the cluster, but much less so in the field: counting stars with 
(g - i) G [0.1, 0.4] and i e [18.5, 19.3] mag (~ 1 Gyr MS- 
TO: 1) and stars with (g-i)o G [0.3, 0.6] and i G [19.4, 20.2] 
mag (~2 Gyr MS-TO: 2), the ratio ATi/JV 2 = 0.21 for the 
field but 0.42 in the direction of the cluster. Hence we iden- 
tify the t ~ 1 Gyr stars with the cluster population, and the 
t ~ 2 Gyr stars with the field population. This is in good 
agreement with the star formation epochs around 1 and 2 
Gyr that Smecker-Hane et al. (2002, their Fig. 5) find for 
another field in the LMC disk. There is also evidence for a 
younger population of several 10 8 yr, both around KMHK 
1603 (Fig. 6, right) and in the disk field studied by Smecker- 
Hane et al. (2002). 

A better match to the optical colour-magnitude dia- 
gram of the cluster stars is obtained by an isochrone of 
slightly younger age, t ~ 800 Myr, yet the best match 
is obtained by an isochrone with t ~ 1 Gyr but with a 
somewhat lower metallicity of [M/H]~ —0.7. Smith et al. 
(2002) also argue that, with [M/H]~ -0.6, LMC red giants 
with Mzams ~ 3 Mq often have metallicities lower than 
what is normally taken as the canonical LMC metallicity, 
[M/H]~ —0.4 to —0.3. Consequently, we also suspect that 
the older, a few Gyr-old, field population has a rather low 
metallicity of [M/H]~ -0.7. 

The near-IR colour-magnitude diagram does not distin- 
guish very well between different intermediate-age popula- 



4 SPECTROSCOPIC EVIDENCE FOR THE 
CARBON-RICH NATURE OF LI-LMC 1813 

Dusty circumstellar envelopes may be separated into carbon- 
rich and oxygen-rich based on their location in IR colour- 
magnitude diagrams, but this diagnostic can be ambigu- 
ous (van Loon et al. 1998; Trams et al. 1999b): the 
colours of LI-LMC 1813, (K - L) ~4.3 mag, (K - 
[12]) = K + 2.51og(Fi 2 /28.3) ~8 mag, and ([12] - [25]) = 
2.51og((F 25 /6.73) x (28.3/Fi 2 )) ~2 mag, are consistent 
with both carbon and oxygen-rich circumstellar material 
(Trams et al. 1999b, their Figs. 7 & 8). To determine the 
dominant chemistry of the stellar atmosphere, we thus re- 
sorted to taking near-IR spectra in the H, K and L-bands. 
These observations are described in full in Appendix A2. 

The L-band spectrum of LI-LMC 1813 (Fig. 9) fea- 
tures strong absorption around 3.1 /Km from HCN and C2H2 
molecules, a narrow absorption feature at 3.55 /im due to 
HCN, and absorption around 3.8 /Km which is attributed to 
C2H2 (Matsuura et al. 2002). These molecules only form in 
great numbers after the supply of oxygen has been exhausted 
in the formation of CO. Hence the strong absorption due to 
carbonaceous molecules is only seen in the spectra of carbon 
stars, that have n(C)/n(0) > 1. The strength of the 3.55 
and 3.8 /jm absorption bands is remarkable, especially as 
they are veiled by continuum emission from the circumstel- 
lar dust at wavelengths A > 3 /Km (see Fig. 12). This hints 
at a contribution from molecules within the dust envelope 
surrounding the star, causing absorption seen against the 
thermal dust continuum emission. 

The carbon-rich nature of LI-LMC 1813 is confirmed 
by the H and K-band spectra (Fig. 10), which clearly show 
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Figure 9. L-band spectrum of LI-LMC 1813. A smooth version 
is offset for clarity. Most of the narrow spikes in the unsmoothed 
spectrum are residuals from the telluric line removal. 
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Figure 11. Equivalent widths of the absorption due to 
HCN+C 2 H 2 at 3.1 pi, W3.1 (circles), and due to HCN at 3.55 
/im, W3.5 (squares), and their ratio, W3.5/W3.1 (triangles), for 
the LMC carbon stars from Matsuura et al. (2002). 
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Figure 10. H and K-band spectra of LI-LMC 1813, on an ar- 
bitrary flux density scale. Smoothed versions of these noisy data 
are plotted offset, and some of the most conspicuous absorption 
features are labelled. 



the absorption bandhead of CO (Av = —2) at 2.293 /tm, 
and the Ballik-Ramsay bandhead of C2 (Av = 0) at 1.77 
fim. The strength of the latter requires a high n(C)/n(0) 
ratio, possibly > 1.4 (Langon &: Wood 2000). A high ratio 
of n(C)/n(0) > 1.2 has also been invoked by Matsuura et 
al. (2002) to explain the high abundance of C2H2 in the 



spectra of LMC carbon stars — including LI-LMC 1813. 
However, the molecular abundances not only depend on the 
free carbon abundance, but also on the local gas temperature 
(Matsuura et al., in preparation). 

Although the equivalent width of the 3.1 fim absorp- 
tion in LI-LMC 1813 is very similar to that in other dust- 
enshrouded carbon stars in the LMC (Fig. 11 and Matsuura 
et al. 2002), the equivalent width of its 3.55 (im absorption is 
the largest of that sample. It also has the reddest (K — [12]) 
colour amongst these stars, which indicates a high mass-loss 
rate. These observations suggest either a highly advanced 
evolutionary state accompanied by a large n(C)/n(0) ratio, 
or a very cool extended atmosphere leading to the abundant 
formation of molecules. 



5 THE MASS-LOSS RATE OF LI-LMC 1813 

5.1 Modelling the Spectral Energy Distribution 

The SED was modelled with the radiative transfer code 
"Dusty" (Ivezic, Nenkova & Elitzur 1999). The input pa- 
rameter "density type" was set to 3, which gives a density 
structure consistent with radiation-driven wind theory. The 
result of the best fit to the data is shown in Fig. 12. 

A stellar temperature of T e g = 2500 K is used as such 
cool photosphere is typical for dust-enshrouded AGB stars 
(van Loon et al. 1998; Groenewegen & Blommaert 1998): 
carbon stars are produced with T c g ~ 2900 K, the superwind 
phase starts at T cS ~ 2700 K, and the coldest carbon stars 
known have T e g ~2200 K (Marigo 2002). Such range in 
temperatures corresponds to a range in derived luminosities 
of only ± 5%, and a range in mass-loss rates of ± 8%. 

The temperature at the inner radius of the dust shell is 
assumed to be Ti nncr = 1000 K as this is the typical tern- 
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Figure 12. Interstellar-extinction corrected spectrophotometric 
energy distribution of LI-LMC 1813, and a Dusty model fit with 
L = 1.5 X 10 4 Lq and ry = 50 (solid). The stellar contribution 
after circumstellar extinction is plotted too (dotted). 



perature at which dust grains condense in a circumstellar 
environment. Since the data used to constrain the SED were 
taken at many epochs, it is quite possible that variations in 
the physical characteristics of the dust shell have occurred 
and therefore modelling with one set of properties will only 
provide an accurate fit for a portion of the SED. Stellar pul- 
sations lead to variability in the luminosity, and therefore 
the dust temperature. A closer fit to the IRAS 25 and 60 
pm data is obtained by assuming Ti nllcr = 800 K. This would 
also increase the stellar contribution in the L'-band, where 
photospheric spectral features are seen. Evidence of a higher 
dust temperature is provided by the much brighter L'-band 
photometric point. Averaged over the pulsational period, a 
value of Ti nllcr = 1000 K seems appropriate. 

As LI-LMC 1813 is a carbon star (Section 4) we use 
amorphous carbon grains (Zubko et al. 1996) ranging in 
size from a = 0.01 to 0.1 /jm according to a power-law 
distribution of the form n(a,)(xa~ q with q = 3.5 (Mathis, 
Rumpl & Nordsieck 1977). We assume the Dusty default 
for the dust grain specific mass density of p g ram = 3 g 
cm -3 . The Dusty model assumes a gas-to-dust mass-ratio 
of f = pgas/pdust = 200. However, intermediate-age stars in 
the LMC may have a higher ratio, as their metallicity can be 
several times lower than solar. This only becomes important 
when calculating the mass-loss rate from the optical depth 
and luminosity. In particular, the radial density distribution 
throughout the radiatively-driven circumstellar outflow de- 
pends on the optical depth but not on the gas-to-dust ratio 
(Elitzur & Ivezic 2001). 

The shell thickness depends on the outflow velocity and 
the length of time for which outflow has occured. The differ- 
ence between the model results adopting r uter /firmer = 10 2 
and 10 4 is slight and only noticeable at A > 30 fxm. As there 
is no reason to assume that the dust shell is either extremely 



thick or thin, we use a relative thickness of 10 3 . The inner 
radius of the dust shell around LI-LMC 1813 is found to 
be runner = 3.9 x 10 12 m, and thus the outer radius lies at 
Pouter = 3.9 x 10 15 m. The Dusty model for the best fit to the 
SED of LI-LMC 1813 predicts an outflow velocity of v = 9.5 
km s , and hence the duration of intense mass loss is es- 
timated at t ~ 10 4 yr. Although the data do not constrain 
this timescale better than by a lower limit of roughly t > 10 3 
yr, it does suggest that LI-LMC 1813 is undergoing a phase 
of prolonged intense mass loss (superwind). 

Dusty produces the shape of the emergent spectrum, 
which is then scaled to match the observed flux. From 
the scaling factor we obtain the bolometric luminosity of 
L = 1.5 x 10 4 L for LI-LMC 1813, or M bo i = -5.72 mag. 
This agrees extremely well both with the original estimate 
of Mboi = —5.78 mag (van Loon et al. 1997), and with the 
value of Mboi,tip-AGB = —5.75 ±0.04 mag as predicted from 
the isochrone fit to KMHK 1603 (section 3.3). 

The optical depth, r, not only determines the shape of 
the SED, but also the mass-loss rate, M, with M oc r. The 
optical depth at visual wavelengths is found to be rv = 50, 
reflecting the highly obscured nature of LI-LMC 1813 (com- 
pare with van Loon et al. 1999b) . The Dusty model returns 
a value for the mass-loss rate valid for L — 10 4 Lq and £ = 
200, with M depending on L and f as MocL 3/4 (£p g r am ) 1/2 . 
At solar metallicity, dust-enshrouded objects have £ ~ 100. 
However, the metallicity of LI-LMC 1813 is ~ 3-5 times less 
than solar (Section 3). As the gas-to-dust ratio scales in- 
versely proportional to the metallicity (van Loon 2000), we 
expect \ ~ 300 to 500 in the dust shell of LI-LMC 1813. 
Hence we estimate the mass-loss rate of LI-LMC 1813 to be 
M = 3.7(±1.2) x 10" 5 M Q yr -1 . 

5.2 Discussion: the superwind of the metal-poor 
carbon star LI-LMC 1813 

LI-LMC 1813 is a carbon star. With its estimated progenitor 
mass of Mzams = 2.2 Mq, 3 rd dredge-up is indeed expected 
to have worked (Marigo 2001) in the absence of HBB (van 
Loon et al. 2001) to raise the photospheric abundance of car- 
bon to above that of oxygen. The observed high n(C)/n(0) 
ratio can be explained in terms of its low metallicity and rel- 
atively high mass as compared to carbon stars in the solar 
neighbourhood that have generally low-mass progenitors of 
Mzams < 2 M Q (see Marigo 2002, her Fig. 6). It also sug- 
gests that LI-LMC 1813 has been a carbon star for already 
quite some time. 

LI-LMC 1813 is indeed a highly evolved AGB star. The 
isochrones that yield the mass estimate (based upon the 
colour-magnitude diagram of the cluster stars) for LI-LMC 
1813 also indicate that its luminosity places it at the very tip 
of its evolution along the AGB. Its luminosity of L — 1.5 x 
10 4 Lq is a few times brighter than the peak in the optically 
bright carbon star luminosity function (Groenewegen 2002), 
but the same as the peak in the dust-enshrouded carbon 
star luminosity function (van Loon et al. 1999b). This is not 
unexpected, as the duration of the TP- AGB is longest for 
stars in the range Mzams ~ 2 to 2.5 M Q (Marigo 2001). 

It is therefore not surprising that LI-LMC 1813 is in 
the superwind phase, experiencing mass loss at a tremen- 
dous rate of M = 3.7 x 10" 5 M yr" 1 . Schroder, Winters 
& Sedlmayr (1999) show that low-mass carbon stars, with 
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Mzams < 1-3 Mq, expell thick shells during a brief period 
of enhanced luminosity as a result of a thermal pulse, but 
that more massive carbon stars such as LI-LMC 1813 suffer 
a longer period of a persistently high mass-loss rate. Our 
analysis of LI-LMC 1813 provides direct evidence for this 
superwind to occur at the tip of the AGB. 

Vassiliadis & Wood (1993) model TP- AGB evolution 
including a prescription for the mass-loss rate based on em- 
pirical data for galactic Mira variables. They argue that the 
maximum mass-loss rate during the superwind phase is set 
by the transfer of momentum from the stellar radiation field 
onto the outflowing matter: M max = v cxp L/c, independent 
of the metallicity or chemical composition. Their prediction 
for the mass-loss rate of M = 1.2-1.4 x 1CT 5 Mq yr -1 during 
the superwind experienced by stars with Mzams = 1.5-2.5 
Mq is three times lower than the actual observed mass-loss 
rate of LI-LMC 1813. However, Gail & Sedlmayr (1986) have 
argued that multiple scattering can enhance the transfer of 
momentum and thus increase the maximum mass-loss rate, 
which was confirmed by observations of dust-enshrouded 
AGB stars in the LMC (van Loon et al. 1999b). 

Wachter et al. (2002) derive a formula for the mass-loss 
rate of solar-metallicity carbon stars, based on dynamical 
models: logM = 8.86 - 1.95 log M/Mq - 6.81 log T cff /K + 
2.47 log L/Lq. According to this formula, the mass-loss rate 
of LI-LMC 1813 is predicted to be M ~ 2+\ x 10~ 5 M Q yr" 1 
for T e ff = 2500 ±300 K. This agrees extremely well with the 
observed mass-loss rate of M = 3.7(±1.2) x 10~ 5 M Q yr -1 . 
Considering that the metallicity of LI-LMC 1813 is three to 
five times lower than solar, the good agreement between the 
actual mass-loss rate of LI-LMC 1813 and the prediction for 
its mass-loss rate at solar metallicity suggests that the mass- 
loss rate — at least during the superwind phase — does not 
strongly depend on metallicity (see van Loon 2000) . 

In her new calculations for AGB evolution and the for- 
mation of carbon stars, Marigo (2002) includes the effects 
of molecule formation within the stellar atmosphere. One of 
the differences between these and previous models is that the 
T e ff drops once molecules are being formed, leading to high 
mass-loss rates of several 10~ 5 Mq yr -1 at n(C)/n(0) > 1.3 
for carbon stars of Mzams > 2 Mq. The observed insensi- 
tivity of M to the initial metallicity may, in the case of 
carbon stars, be due to enhanced formation of carbonaceous 
molecules as a result of a higher n(C)/n(0) ratio at lower 
metallicity. Alternatively, metallicity insensitivity of the pul- 
sation mechanism has been suggested as a possible reason 
for the metallicity insensitivity of the mass-loss rate during 
the superwind phase (van Loon 2002). 

When an AGB star arrives on the Mira period- 
luminosity (P-L) sequence, its mass-loss rate has become 
so severe that further evolution is along tracks of constant 
luminosity, with the radius (and hence pulsation period) 
growing as the structure of the mantle re-adjusts itself to 
the decrease in mass (e.g. Wood 1998). Describing the Mira 
as a polytrope, RtyM = constant, and the pulsation as a 
harmonic oscillator, 2n/P ~ l/tf re efaii — g(R)/R, the period 
evolves as Pocl/M (van Loon 2002). 

If LI-LMC 1813 has evolved from P ~400 d on the Mira 
P-L (Mzams = 2.2 Mq) to P ~680 d at present, it has a 
current mass of only M = 1.3 Mq. According to Marigo 
(2001), its final mass will become M W d = 0.6 M Q . This 
means that, although it has already lost AM ~ 0.9 Mq, it 



must still lose another AM ~ 0.7 Mq. Stars of Mzams = 2.2 
M with correspondingly long periods (P ~ 1500 d) have 
not (yet) been found, the longest known period for such star 
being P ~900 d (Wood 1998). This is not surprising, as it 
would perhaps take as little as 10 4 yr to lose the remainder 
of the stellar mantle, rendering these objects extremely rare. 



6 SUMMARY 

As a result of a large observational effort, the mid-IR point 
source LI-LMC 1813 is now one of the few AGB stars expe- 
riencing the superwind phase of intense mass loss for which 
the progenitor mass, initial metallicity, atmospheric chem- 
istry, luminosity and mass-loss rate are accurately known. It 
thus provides an important testcase against which to gauge 
models for the evolution and mass-loss of AGB stars. 

Absorption features due to carbonaceous molecules in 
the near-IR spectra of LI-LMC 1813 indicate that it is a 
carbon star. From its location in the LMC cluster KMHK 
1603, its progenitor mass is inferred to be Mzams = 2.2 Mq, 
which is typical for a carbon star, whilst its initial metallicity 
is three to five times lower than solar. 

The object is heavily obscured by its circumstellar dust, 
which shines brightly at IR wavelengths. From the modelling 
of the spectral energy distribution, a mass-loss rate of M = 
3.7 x 10~ 5 Mq yr -1 is estimated. This is very high, especially 
for such metal-poor star. As the luminosity of LI-LMC 1813 
(L = 1.5 x 10 4 Lq) places it at the very tip of its AGB, 
this observation provides direct evidence for the superwind 
to represent the final episode in the AGB evolution. 

Although the observed mass-loss rate for LI-LMC 1813 
is a few times higher than that predicted by the formalism in 
Vassiliadis & Wood (1993), it is in very good agreement with 
the prescription based on the dynamical models of Wachter 
et al. (2002) for solar metallicity. This suggests that the 
mass-loss rate in the superwind phase may be metallicity in- 
sensitive (van Loon 2000, 2002), perhaps because of a larger 
abundance of carbonaceous molecules other than carbon- 
monoxide in metal-poor stellar atmospheres (Matsuura et 
al. 2002; Marigo 2002). 
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APPENDIX A: OBSERVATIONS 
Al Photometry 

Al.l Optical imaging (Dutch 0.9m, La Silla) 

The direct imaging camera at the Dutch 0.9m telescope at 
La Silla, Chile, was used on the six nights of December 25 
to 30, 1996, to obtain deep images of a 3.77' x 3.77' region 
around the cluster KMHK 1603, through Gunn g (A = 5148 
A, AA = 81 A), r (A = 6696 A, AA = 103 A) and i 
(A = 7972 A, AA = 141 A) filters (Thuan & Gunn 1976; 
Wade et al. 1979). The total integration time amounts to 4 
hours per filter, split into 5-minute exposures to avoid satu- 
ration and to allow for refocussing in order to reach the best 
image quality. The pixels measure 0.442" x 0.442" on the 
sky, and stellar images on the combined (shift-added) frames 
have a Full Width at Half Maximum (FWHM) of ~ 1.5", 
though some individual frames show stellar images of ^1.1" 
FWHM. The CCD frames were reduced using standard pro- 
cedures within the ESO-MIDAS package. The reader is re- 
ferred to van Loon et al. (2001) for more details. 



2MASS and Denis is transformed onto the SAAO system 
(Carter 1990; Carpenter 2001) before further analysis. 



A1.4 L' -band imaging (ISAAC, Paraxial) 

The Infrared Spectrometer And Array Camera (ISAAC) at 
the ESO VLT at Paranal, Chile, was used on October 19, 
2000, to image LI-LMC 1813 through an L'-band filter (A = 
3.78 (im, A A = 0.58 /jm). Because of the high background, 
the image was obtained after chopping and nodding (both by 
10" in Declination), with 9 sub-integrations of 0.104 second 
at each chop position, and a total on-source integration time 
of 1 minute. The field-of-view of 72" x 72" (at 0.071" pixel" 1 ) 
ensured that KMHK 1603-1 was imaged simultaneously with 
LI-LMC 1813. Conditions were very good, with 6% relative 
humidity at a temperature of 15° C. The seeing was 0.67" at 
Zenith, with the target at an airmass of 1.442. 

The data were reduced using ESO-MIDAS and the 
Eclipse package (Devillard 1997). First the electrical ghosts 
were removed (lS_GHOST in Eclipse) . Then the relative pixel 
responses were calibrated using a fiatfield image, and aper- 
ture photometry was performed on the difference of the two 
images at each nodding position. The photometry was cali- 
brated against that of the G3-giant standard star HR 2354 
ill = 4.989 mag; van der Bliek, Manfroid & Bouchet 1996: 
virtually identical filter as that of ISAAC) at similar airmass 
(1.392) and only 10 minutes after LI-LMC 1813. 



A1.2 J and K-band imaging (IRAC-2, La Silla) 

The near-IR camera IRAC-2b at the ESO/MPI 2.2m tele- 
scope at La Silla, Chile, was used on January 2, 1996, to 
image the field around the IRAS point source LI-LMC 1813. 
Lens C was chosen to cover a field of view of 133" x 133", 
with a pixel scale of 0.49". A sequence of 12 images in the K- 
band filter were obtained, each consisting of ten 3-second ex- 
posures, shifted by 5" in Right Ascension (RA) with respect 
to the previous image. This procedure was then repeated 
backwards, using the J-band filter. The image is deepest in 
the central ~ 75" x 130". The photometry is transformed 
onto the SAAO system (Carter 1990). The reader is referred 
to van Loon et al. (1997) for more details. 

A 1.3 Near-IR imaging (2MASS & DENIS) 

Archival J, H and K-band data is available from the sec- 
ond incremental data release of 2MASS. It contains 8 stars 
within 1' of LI-LMC 1813 which were detected in the K 
band down to A'nmit ~ 15 mag, including LI-LMC 1813 it- 
self, on 26 October 1998. Especially valuable is the detection 
of LI-LMC 1813 in the H band, which is thereby the short- 
est wavelength at which the object has been seen. The Denis 
point source catalogue towards the Magellanic Clouds (Cioni 
et al. 2000) provides additional I, J and K s -band data, but 
its sensitivity only reaches K Bi ij m j t ~ 12 mag and no star 
was detected within 1' of LI-LMC 1813 except for the opti- 
cally brightest star in the cluster, KMHK 1603-1, which was 
detected in the I and J-bands. The photometry from both 



A 1.5 Mid-IR imaging (MSX) 

LI-LMC 1813 is included in the Mid-course Space experi- 
ment (MSX) Point Source Catalogue (Price et al. 2001). It 
was only detected in band A (effective wavelength Ao = 8.28 
fim). The transmission of this band runs from A ~ 6.1 to 10.8 
/jm, though, peaking at A pca k ~ 9.9 fim, and the measured 
flux density thus depends heavily on the spectral slope of the 
object — which is greatly affected by circumstellar emission. 



A1.6 N-band imaging (TIMMI-2, La Silla) 

The TIMMI-2 instrument at the ESO 3.6m telescope at La 
Silla, Chile, was used on January 19, 2001, to obtain an 
image of LI-LMC 1813 through the N2 filter (A = 10.6 /nm, 
with a transmission > 50% between A = 9.74 and 11.33 /im). 
At a pixel scale of 0.3" pixel" 1 , the slightly vignetted field 
measures 96" (RA) x 72" (Dec). The standard technique for 
observing in the thermal IR was employed, chopping 20" in 
Dec and nodding 20" in RA. The total on-source integration 
time was 15 minutes. The images were combined, and the 
photometry was obtained, using standard commands within 
ESO-MIDAS. The K3 giant a Hya was used as a prime 
calibrator, adopting a flux density of Fn2 = 129 Jy (Cohen 
1998). The photometry was checked against the IRAS 12 
/jm flux densities for the K2 and Kl giants A2 Tuc (F12 = 3 
Jy) and 7 Pic (F 12 = 7 Jy). 
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A 1.7 Mid-IR photometry (IRAS) 

Data at 12, 25 and 60 /mi were retrieved from the IRAS 
data base server in Groningen 1 (Assendorp et al. 1995). 
The Gipsy data analysis software was used to measure the 
flux density from a trace through the position of the star 
(Gipsy command SCANAm). For a discussion about the sen- 
sitivity and reliability of this procedure for point sources in 
the LMC, see van Loon et al. (1999b). 

A2 Spectroscopy 

A2.1 Optical spectroscopy (EFOSC-2, La Silla) 

The EFOSC-2 instrument at the ESO 3.6m telescope at La 
Silla, Chile, was used on October 31, 2002, to obtain an op- 
tical spectrum of KMHK 1603-1 between A = 4700 and 6700 
A with grism #9. The slit was narrowed to a width of only 
0.5" projected on the sky, in order to boost the spectral re- 
solving power to R ~ 2000. Two spectra were taken, with ex- 
posure times of 300 seconds each, which were then combined. 
The bias level was subtracted first. Differences in illumina- 
tion and responsivity across the image were corrected using 
flatfield spectra of an internal quartz lamp. The wavelength 
axis was calibrated using two exposures of a He/Ar lamp, 
with an accuracy to well within the spectral resolution. The 
sky-subtracted spectra were corrected for extinction by the 
Earth's atmosphere, and for the response function as derived 
from a spectrum of the HST standard star HD 49798. 

A2.2 H and K-band spectroscopy (OSIRIS, CTIO) 

The Ohio State InfraRed Imager/Spectrometer (OSIRIS) at 
the 4m Blanco telescope at Cerro Tololo (CTIO) , Chile, was 
used on 6 July, 2002, to obtain a spectrum of LI-LMC 1813 
covering the 0.9 to 2.4 /jm wavelength range. The instrument 
was used first in imaging mode to acquire the star in the slit. 
The spectra were acquired in cross-dispersed mode through 
a 1" x 30" slit, providing a resolving power of R = 1200. Four 
spectra were taken, with each time the star moved over 5" 
along the slit. The slit orientation was East- West. 

The data reduction was performed using Iraf. In the 
individual spectra bad pixels were eliminated and flatfield 
and linearisation corrections were applied. Since the spec- 
tra were taken in the morning almost at sunrise, the sky 
brightness changed considerably during the course of the 
observations. From each image another image (the nearest 
in time) was subtracted to perform a first sky subtraction. 
The individual spectra were then extracted and combined. 
Stars of spectral type A and B were observed as close to 
the LI-LMC 1813 spectrum's airmass (1.756) as possible, 
in order to correct for telluric absorption features and the 
instrumental response. The raw spectra of these reference 
stars were obtained in the same manner as LI-LMC 1813. 

The IRAS data base server of the Space Research Organisation 
of the Netherlands (SRON) and the Dutch Expertise Centre for 
Astronomical Data Processing is funded by the Netherlands Or- 
ganisation for Scientific Research (NWO). The IRAS data base 
server project was also partly funded through the Air Force Of- 
fice of Scientific Research, grants AFOSR 86-0140 and AFOSR 
89-0320. 



The hydrogen and helium series appearing in these spectra 
were removed by linear interpolation. The cleaned spectra 
of the early-type stars were divided by a blackbody curve 
to eliminate the slope due to their intrinsic energy distribu- 
tions (see Messineo et al., in preparation, for more details). 
Temperatures were adopted from Gray (1992) and Vacca, 
Garmany & Schull (1996). The wavelength calibration was 
obtained using the telluric OH lines (Iraf tasks refspectra 
and DISPCOR; see Oliva & Origlia 1992). The spectra were 
dereddened (Iraf task deredden) using the extinction law 
of Cardelli, Clayton & Mathis (1989) — essentially the same 
as that of Mathis (1990). 

A2.3 L-band spectroscopy (ISAAC, Paranal) 

ISAAC at the ESO VLT at Paranal, Chile, was used on 
December 14, 2001, to take an L-band spectrum of LI-LMC 
1813. The total on-source exposure time was 36 minutes. 
With a slit width of 0.6", the resolving power was R ~ 600. 
The optical seeing was 1.2", but the infrared seeing was 
~ 0.5". The background was subtracted by chopping with 
a chopping throw of 15". Telescope nodding and jitter was 
used to correct for bad pixels. The data was reduced using 
the Eclipse package (Devillard 1997). Variations of the pixel 
response across the array were corrected with a flatfield of 
the twilight sky. The wavelength calibration is based on the 
exposure by an Ar+Xe arc lamp with the same wavelength 
setting as the target observation. 

The weather was photometric. The B9 IV star HIP 
029729 was used as a telluric standard and flux calibration 
star. For the flux calibration we use the effective temper- 
ature, T e ff, and (K — L) colour for a spectral type B9 V 
instead of B9 IV, because no accurate value for T e g is avail- 
able for sub-giants. The spectrum is represented by a black 
body with T cS = 10,700 K (Tokunaga 2000). The K-band 
magnitude of this star is K — 7.536 mag (van der Bliek et 
al. 1996) and (K - L) = -0.03 mag for a B9 V star, hence 
the L-band magnitude is estimated to be L — 7.566 mag. 
The L-band magnitude is converted into Jansky according 
to the zero- magnitude in van der Bliek et al. (1996). 



APPENDIX B: THE OPTICALLY BRIGHT 
STAR KMHK 1603-1 

The optical and infrared colours and magnitudes of the 
bright star KMHK 1603-1 are consistent with those of an 
F-type supergiant in the LMC. It could then either be a 
field supergiant of relatively young age (t <C 10 8 yr), or a 
post-AGB star in the cluster KMHK 1603. The latter case 
would be especially interesting, as we could then compare 
the properties of an AGB and a post-AGB star of the same 
main-sequence mass and metallicity. However, it could also 
be a galactic foreground main-sequence star. 

Modelling of the photometry (Fig. Bl) suggested a tem- 
perature of (roughly) T e s ~ 6500-7000 K. If at the distance 
of the LMC, it will then have a bolometric luminosity of 
L ~ 1.5 x 10 4 L Q , i.e. virtually identical to that of LI-LMC 
1813. As the luminosity of a post-AGB star (before reaching 
the White Dwarf cooling track) is the same as it had when it 
left the AGB, this lended further support for the hypothesis 
that KMHK 1603-1 is a post-AGB star in the KMHK 1603 
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Figure Bl. Spectro-photometric energy distribution of KMHK 
1603-1, corrected for interstellar extinction, assuming that it is 
located in the LMC. Blackbodies of T = 6500, 7000 and 7500 K 
are shown for comparison. 
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Figure B2. Ha and H/3 line profiles in the optical spectrum of 
KMHK 1603-1. Peaking to within ~30 km s" 1 from the helio- 
centric rest velocity, they rule out LMC membership. 

cluster. Therefore, an optical spectrum was taken in order 
to confirm or reject this. 

The Ha and H/3 photospheric absorption lines in the 
spectrum of KMHK 1603-1 peak within ~30 km s~ from 
the heliocentric rest velocity (Fig. B2), ruling out LMC 
membership. On the basis of the relative strengths of the 
Ha and H/3 lines, A5174 Mg I triplet (Pritchet & van den 
Bergh 1977), AA5890, 5896 Na I D doublet and the A5270 



and A6495 Fe I lines, and by comparison with the spectral 
library of Jacoby, Hunter & Christian (1984), the spectral 
type is estimated to be F6 V (T e ff = 6400 K). This suggests 
a bolometric luminosity of only L ~ 15 L.0, and hence a dif- 
ference in the distance modulus of A(m — M) ~ 10 mag, 
which corresponds to a distance of d ~ 500 pc, or a distance 
above the galactic plane of z ~ 250 pc. From a comparison 
of the acquisition image of October 2002 with the imaging 
in December 1995, no proper motion could be detected in 
excess of /i ~ 0.02" yr _1 , but the proper motion for a star 
at d = 500 pc and galactic longitude of in = 277° is also not 
expected to exceed this upper limit. 
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